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Abstract 

By using quark Yukawa matrices only, we can construct renormalization invariants that 
are exact at the one-loop level in the standard model. One of them G is accidentally consis¬ 
tent with unity, even though quark masses are strongly hierarchical. We calculate a lepton 
version of the invariant P for Dirac and Majorana neutrino cases and hnd that P can also 
be close to unity. For the Dirac neutrino and inverted hierarchy case, if the lightest neutrino 
mass is 3.0 meV to 8.8 meV, an equality U = P can be satished. These invariants are 
not changed even if new particles couple to the standard model particles, as long as those 
couplings are generation independent. 


1 Introduction 

In the standard model (SM) of particle physics, fermion masses are determined by their Yukawa 
couplings to the Higgs field. These couplings are free parameters and cannot be predicted within 
the SM. There may be some principles that determine the couplings at a high energy scale, such 
as a grand unified theory or a flavor symmetry. To compare the predictions that come from those 
high energy theories and observed values by experiments, we have to consider the renormalization 
group evolution of the parameters. Renormalization group invariants are useful because we can 
infer the high energy physics with low energy inputs. Invariants that are constructed only by 
quark Yukawa coupling matrices were reported in Ref. [1]. One of them, which we call J'^, is 
close to unity. This fact seems surprising because their building blocks are far from one. This 
invariant is not changed even if there are some new particles that contribute to the beta functions 
in a generation-independent way since those contributions are canceled (as in Ref. P; see the 
Appendix). For example, if right-handed neutrinos exist and have Yukawa couplings of 0(1), they 
largely change the evolution of other Yukawa couplings [2], but do not affect J'^. 

Today, we know that neutrinos have masses, so it is likely that they are coupled to the Higgs 
field as other fermions. There are still many mysteries in neutrinos: their type (Dirac or Majorana), 
the mass hierarchy (normal or inverted), and the absolute masses. We calculate a renormalization 
group invariant P in the lepton sector that is similar to P to see implications on these mysteries. If 
neutrinos are Dirac type, the relation between the neutrino Yukawa coupling matrix and their mass 
matrix is the same as other fermions. If neutrinos are Majorana type, we need some assumption to 
relate these matrices. In this paper, we consider simple cases of the type-/ seesaw mechanism [3H5]. 
To calculate P, we introduce a useful parametrization for neutrino Yukawa couplings. 

As a result, we hnd that P can also be 0(1) regardless of the type of neutrinos. For the 
Dirac neutrino and the inverted hierarchy case, an equality between quark and lepton sector 
jq — ji hold if the lightest neutrino mass is 3.0 meV< mnghtest < 8.8 meV. For the Majorana 
neutrino cases, the invariant depends on many unknown parameters. Not only P = P, but also 
an interesting condition P = 1 can be realized for Majorana cases. 
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2 Renormalization group invariants 


The Yukawa terms of the Lagrangian that we consider in this paper are 


-L = + IrYi^^Lj^ + H.c.. 


( 1 ) 


We assume there are three right-handed neutrinos Renormalization group equations (RGBs) of 
Yukawa couplings mu (see the Appendix) can be simplihed by rewriting with Hermitian matrices 
«/ = Pn. (/ = u, d, u, 1). These quantities do not depend on the bases of right-handed particles 
Jr. There is a combination of the quark Yukawa matrices that is exactly invariant under the 
renormalization group evolution at the one-loop level [1], 


/'' = (Det[a„arf]) 


Tr[a„ad] 

Tr[(a„a;d)-i]‘ 


( 2 ) 


This invariant does not depend on the basis of particles in the flavor space, since a unitary 
transformation of the basis of the left-handed quark doublet Ql is also canceled by taking traces 
and determinants. is invariant even if we introduce new particles, as long as they couple to the 
SM particle in a generation-independent way (see the Appendix). 

In the basis that is diagonal, we can express the Yukawa matrices as 


= Diag[m^, ml, ml]/v^, (3) 

ad = 1/Diag[m^, ml, ml]Vyv^. (4) 


V denotes the quark mixing matrix and v is the vacuum expectation value of the Higgs doublet 
<h. By substituting the observed values (renormalized to Mz) [Hl[9], three factors in Eq. ([2]) are 

Det[a„ad] = (jriumcmtmdmsmbY= 6.26 x 10“^®, 

Tr[a„ad] = ^ Vgi| = 2.63 x 10"^, 

Tr[(a„ad)"^] = ^ V = 5.76 x 10^®. 

The summations are over jS = u,c,t and i = d,s, b. These factors are much deviated from unity. 
The invariant is, however, close to unity, 

p = 1.35l|j = 0.72 ... 1.89 (5) 


[Icr errors are used in this paper except Eq. (|9])]. The large errors of P come from the errors 
of quark masses < 30% [8]. It is surprising that this invariant is so close to and consistent with 
unity. By keeping leading terms and an approximation V ~ 1, ~ 1 means 


mumtmdmb 

- ~ 1 . 

mlm? 


( 6 ) 


Such a mass relation was pointed out in Ref. m. In this paper, we calculate the invariant in the 
lepton sector. 


f = 


(Det[aj.ai]) 


Tr[a^a;] 
Tr[(aj.az)-i] ‘ 


( 7 ) 


We can see that P is invariant in a manner discussed in Ref. [T]. The relation between neutrino 
Yukawa couplings and the observed neutrino mass matrix depends on the type of neutrinos: Dirac 
or Majorana. 
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3 Dirac neutrinos 


If neutrinos are Dirac type, the neutrino Yukawa matrix is directly converted to neutrino mass 
matrix rrii, as 

= -m^ = -URmU\ rriij = mAj {i,j = 1,2,3) 

V V 

U and Ur are unitary matrices that diagonalize rrii,. In the basis where the charged lepton mass 
matrix is diagonal, U corresponds to the neutrino mixing matrix. The invariant P can be explicitly 
written similar to U, 


P 


, ^-4/3 

TUi \Uj3i\^ 


( 8 ) 


The summations are over /3 = e, p, r and i = 1, 2, 3. Among the parameters in Eq. ([ 8 ]), absolute 
values of the neutrino masses cannot be determined by oscillation experiments. The strongest limit 
on the neutrino mass scale is set by the observation of the cosmic microwave background 


< 0.23 eV (95% limit) 


(9) 


so we consider the region < 0.1 eV in this paper. The dependence of P on the lightest neutrino 
mass is shown in Fig. [TJ We have used charged lepton masses (renormalized to Mz [S]) and the 
results of neutrino oscillation experiments [12] • To better understand Eq. ([H|), we discuss three 
limiting cases below. 


(i) Degenerate masses 

The simplest case is that neutrino masses are almost the same. If mi = 0.1 eV, by using 
mass differences given in Ref. [12], we obtain m 2 = 0.1004 eV and m 3 = 0.11 eV for the 
normal hierarchy (NH) case and m 3 = 0.087 eV for the inverted hierarchy (IH) case. By 
approximating mi ~ m 2 — m 3 and noting that the charged lepton masses are hierarchical 
me -C m^ -C mr, we can roughly estimate 


tI ( 3 1 - 4/3 rnjmlJ2^\Uri\^ 

I ~ (mimem^m,) ' ^ 

mi me^2^j If/eil 


> 2/3 

m^m^ \' pi 

I— -'dJ 


mf 


= 0.18648 ±0.00001. 


( 10 ) 


Here we have used the unitarity of U. Interestingly, P approximately does not depend on 
the neutrino masses, and this feature can be seen in the region of mi ~ 0.1 eV of Fig. m 


(ii) Normal Hierarchy (NH) 

We take mi m 2 < m 3 . By subtracting the most dominant term of P, 

^ f rriemr mimg ^ |f/r3p _ 

\ m^ m\ ) |f4iP '^\m\ 

We have used the conventional parametrization of U given in Ref. [12]. The experimental 
error of P is roughly estimated from this equation to be 195 %. 


2/3 9 


-23 


-13 


( 11 ) 
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^lightest/ 


Figure 1: Neutrino mass dependence of the renormalization group invariant P for the Dirac 
neutrino case. Upper and lower curves show the IH and NH cases with la deviation roughly 
estimated from Eqs. flTTl) and flT 2 l) . The straight line shows the lower limit of I^. 


(iii) Inverted Hierarchy (IH) 

We take m 2 S> m 3 . In this limit, m 2 = 0.0495 eV and mi = 0.0487 eV, so we can approxi¬ 
mate mi ~ m 2 . The dominant term of is 


/ |UrlP + \Ut2? ^ tI f ^ 

V m2 m2y lUesP “ s?3 


( 12 ) 


Since S 13 <C 1, P tends to be larger for the IH case than that for the NH case. The 
experimental error of P is roughly estimated from this equation to be ^ 31 %. 


From Fig. [H we hnd that P is 0(0.1) 0(1) for both mass hierarchy cases. In the IH 

case, P = can be satished for narrow and hierarchical regions of the lightest neutrino mass 3.0 
meV< m 3 < 8.8 meV mi,m 2 (see Tabled]). This mass range gives 

0.101 eV < ^ mi < 0.109 eV. (13) 

Such a total neutrino mass is within the reach of future cosmological and astrophysical surveys [13]. 
Unlike J'^, an interesting condition P = 1 is highly excluded for both NH and IH cases. 


4 Majorana neutrinos 

If neutrinos are Majorana type, we need some assumption to relate their mass and the Yukawa 
matrix. Here we assume that the Majorana masses are generated through the type-I seesaw mech- 
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anism [3HS]. We add a Majorana mass term of three right-handed neutrinos to the Lagrangian, 

-I- , Mjj = MjSij {I, J = 1, 2, 3). (14) 

We take the basis of ur such that M is diagonal. At a renormalization scale larger than the 
heaviest right-handed neutrino, RGBs of Yukawa couplings are the same as the Dirac neutrino 
case [13], so P is invariant in this region. We expect that the Majorana masses of ur are much 
larger than the Dirac mass term. Then we obtain the neutrino mass matrix 

(15) 


In this case, we can parametrize the neutrino Yukawa couplings as [T5| 

= U^/mR^ MRy/mUyv"^. (16) 

R is an arbitrary complex matrix that satisfies R^R = 1. This matrix cannot be determined by 
low energy experiments. The determinant of Eq. (ll 6 p does not depend on R, 

Det[a^] = I)et[mR^MR/v^] 

= mim2ra'iMiM2M‘i/. (17) 


Here we have used Det[i?l]Det[i?]= (±1)^ = 1. In the following, we assume degenerate masses 
Ml = M 2 = M 3 = Mtv. By these parametrization, the invariant is 


, , ^ ^ Tr rRlR-y/mt/lmff/Yml 


(18) 


P does not depend on M explicitly, but depend on R^R. If R is real, R^R = R^R = 1. The 
unitary matrix U have freedom of Majorana phases Ai^s, 


U = U^Dm, Dm = Diag[e*"b 1, 


iXi 




(19) 


where is the mixing matrix same as the Dirac neutrino case. Note that Dm and m are 
commutative, so if R is real. 


DM\/m\/mD^j^U^yv‘^ = U^mU^yv'^, 


( 20 ) 


which means that P does not depend on the Majorana phases. 

The matrix R^R in Eq. fflSj) is Hermitian and orthogonal. We introduce a useful parametriza¬ 
tion for such a matrixji] 


R^R = VlT^I - - 1 )^ + lA, 


0 CI3 — 0,2 

A = I —03 0 Oi 

02 —Oi 0 


, d~ = (ai, 02, ^3), o = |o| 


( 21 ) 

( 22 ) 


^ Another possible parametrization is — exp(iI3), where B is a real antisymmetric matrix |16| . This is 
concise, but values of each component are not explicit. 
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If R is real, R^R = R^R = 1, which corresponds to a = 0@ In principle, a, can take any real 
value, and P can be inhnitely large. Consider the case 02 = 03 = 0, Oi = a 3> 1. In this limit, we 
obtain 


/O 0 0\ 

R^Rc^a 0 1 i . (24) 

Vo 1/ 

Trace parts in Eq. flTSll do not depend on mi, so the invariant is 

P oc (25) 

P diverges in the limit mi —)■ 0 . The limit a ^ 1 is, however, highly unnatural, since it means that 
the absolute values of the components of are much larger than that of YjY^y [compare Eq. 
(IT^ and ffTOj) ]. To be concrete, we consider the region |aj| < 5. Note that a small imaginary part 
of R is enough to produce the baryon asymmetry of the Universe through leptogenesis [TTlfTS] . 

The dependence of P on the lightest neutrino mass is shown in Figs. [2]and[3l For each NH and 
IH case, we have shown three curves: (i) a* = 0, which corresponds to the case of real R, (ii) /^(O.l 
eV) is largest and (hi) smallest. Two Majorana phases are set to zero in the (ii) and (iii) cases. We 
have used charged lepton masses renormalized to Mz [ 8 ] and neutrino oscillation parameters [ 12 ] . 
P is dehned at the energy scale larger than Mjv, but we can expect that the renormalization 
group effect between Mz and M^v is negligible. The reasons are as follows. First, P include 
charged lepton masses in the form of ratios, which are almost constant in renormalization group 
evolution. Second, P does not depend on the evolution of the overall scale of neutrino masses, 
and the renormalization effect on mixing angles is negligible [T9II2T] since the Yukawa coupling of 
the T lepton is small in the SM (the effect is estimated to be smaller than 1° [20]). 

As we can see from Figs. [2]and[3l the invariant P largely depends on the imaginary part of 
R^R. If R is real, the equality P = R is satished only in the IH case with small m 3 . 


0.19 meV < m 3 < 1.5 meV 


( 26 ) 


(see Tabled]). Majorana neutrino masses can be searched by not only cosmological observations 
but also neutrinoless double beta decay experiments that are sensitive to m^^ = | 

Eq. we obtain 


0.098 eV <^mi < 0.100 eV, (27) 

0.018 eV <mpf} < 0.048 eV (28) 


(Majorana phases are considered for mpg). These regions can be searched in the foreseeable 
future [T3l|22]. If R is complex, P can be much larger, and an interesting condition P = 1 can be 
satished for both NH and IH cases. 

^Without this condition, a general orthogonal Hermitian matrix can be written as 


± 


\/l + a^I — {p 1 + ± l)dd^/a^ + iA 


(23) 
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Figure 2: Majorana neutrino, NH case. The blue solid curve is the case of real R, or a = 0. 
The upper and the lower dashed curves are cases in which /*(0.1 eV) is largest and smallest 
[a^ = (5.00, 2.89, 0.572) and = (—1.98, 3.80, —5.00)]. The lower limit of the invariant R is also 
shown by a straight line. 



5 Summary 

We can construct renormalization group invariants with Yukawa coupling matrices only. As 
pointed out in Ref. [T], an invariant R in the quark sector is close to one, even though quark 
masses are highly hierarchical. If this fact implies a rule for the structure of Yukawa matrices at 
the high energy scale, it would also hold in the lepton sector. In this paper, we have calculated 
an invariant P of the lepton Yukawa couplings, which is dehned similar to R. Note that R and 
P remain constant even if new particles couple to the SM particles generation independently. 
To calculate P for Majorana neutrino cases, we introduced a useful parametrization of neutrino 
Yukawa couplings. We found that P can also be 0(1) for both neutrino types. As the quark sector 


Table 1: Summary of our numerical result. For Majorana cases, R is assumed to be real, /peak 
denote the extremum of P (when the lightest neutrino mass is mughtest = ’^peak)- mp^jq is the 
lightest neutrino mass that satishes the condition P = R within one standard deviation. 



Dirac 

Maj 

orana 


NH 

IH 

NH 

IH 

■^peak 

U.ZUU_Q Q25 

u. /ou_g oei 

n iqi+u.uiu 

u.iai_o.oi8 

U. 061 

mpeak/meV 

7.06 

5.30 

7.59 

0.568 

mp^iq /meV 

- 

3.0... 8.8 

- 

0.19... 1.5 
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Figure 3: Majorana neutrino, IH case. The upper and the lower dashed curves are taking = 
(5.00, 2.60, —0.235) and = (—3.58, 5.00, —5.00). Straight horizontal lines show the central value 
and lower limit of . Other notations are the same as in Fig. [2l 



invariant, P = 1 is consistent for Majorana neutrino cases, but not for Dirac neutrino cases. A 
hypothetical equality P = /'^ tends to favor inverted and strongly hierarchical neutrino masses, 
and such a mass spectrum can be explored by future experiments. 
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Appendix Derivation of invariants 


In this Appendix, we see that /'^ is invariant in the SM and some extensions of it. The invariance 
of P can be checked in a similar way. The RGBs of the quark Yukawa couplings are 


d 3 

2 ^ulNgXNgl 

d 3 

^ IdlNgXNg- 


(29) 

(30) 


dt = (i/i/(167r^/x), where /x is the renormalization scale. /wgxWg is a unit matrix of the size Ng{= 3 
in the SM) and 'ju^ are flavor-independent terms of the beta functions. They come from the wave 
function renormalization of the Higgs doublet <F and radiative corrections by gauge interactions. 









The RGBs of basis-independent quantities are calculated as 


^Tr[a„Q!rf] = 2(7„ + 'jdf^^iauad], 

(31) 

^Tr[(a„ad)"^] = -2(7^ J-7rf)Tr[(Q!„Q!d)“^], 

(32) 

^Det[a„ad] = + 7d)Det[a„ad]. 

(33) 


Then, the running of is 

^ (Det[a„ad]“^/^»Tr[a„arf]Tr[(Q!„arf)"^]"^) 

= (-^2iV, + 2-(-2)^ hu + 7<i)I^ 

= 0, (34) 

so is a renormalization group invariant. We can easily see that Tr[Q!„arf]Tr[(a„ad)“^] is also 
invariant. Our derivation does not depend on the content of 'yu,d- R means that new particles that 
couple only to gauge fields or Higgs doublet do not change /^. If new particles couple to quarks 
generation-independently, their contribution are also absorbed in 7^ rf, so remain invariant. For 
an example, right-handed neutrinos do not couple to quarks but couple to the Higgs doublet. 
They add TT[a^]lNgxNg to the right-hand sides of Eqs. fl2^ and fl30l) . Such terms can largely 
change the evolution of W and Yd |2], but P remain invariant. 

The essential reason that we could find those invariants is that the generation-dependent 
terms can be canceled, Y~^-^Yu + Y^^-^Yd oc IwgxNg- Two Higgs doublet models and minimal 
supersymmetric models change the Yukawa sector of the SM (see Ref. [7] for RGBs), and then we 
cannot construct invariants in a parallel way. 
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